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We explore some particle physics implications of the growing evidence for a helical primordial 
magnetic held (PMF). From the interactions of magnetic monopoles and the PMF, we derive an 
upper bound on the monopole number density, n(to) < 1 x 10~^°cm“®, which is a “primordial” 
analog of the Parker bound for the survival of galactic magnetic helds. Our bound is weaker 
than existing constraints, but it is derived under independent assumptions. We also show how 
improved measurements of the PMF at different redshifts can lead to further constraints on magnetic 
monopoles. Axions interact with the PMF due to the ga~tpE ■ B/47r interaction. Including the 
effects of the cosmological plasma, we Hnd that the helicity of the PMF is a source for the axion 
held. Although the magnitude of the source is small for the PMF, it could potentially be of interest 
in astrophysical environments. Earlier derived constraints from the resonant conversion of cosmic 
microwave background (CMB) photons into axions lead to ga-^ < 10“® GeV“^ for the suggested 
PMF strength ~ 10“^^ G and coherence length ~ 10 Mpc. Finally we apply constraints on the 
neutrino magnetic dipole moment that arise from requiring successful big bang nucleosynthesis in 
the presence of a PMF and we hnd 


I. INTRODUCTION 

There is growing evidence for the existence of an in- 
tergalactic magnetic held from the observation of high 
energy gamma rays. It is likely that a magnetic held in 
intergalactic space would have been created in the early 
universe, since astrophysics alone is not expected to gen¬ 
erate helds on such large length scales. (For a recent re¬ 
view on cosmic magnetic helds, see Ref. [1].) The discov¬ 
ery of a primordial magnetic held (PMF) has important 
ramihcations for cosmology as it allows one to test mod¬ 
els of magnetogenesis, which are often tied to the physics 
of inhation [2] cosmological phase transitions [3, 4], and 
baryogenesis [5-7]. The presence of a PMF after cos¬ 
mological recombination can also aid in the formation of 
hrst stars [8] and provide the seed held for the galac¬ 
tic dynamo [9]. Additionally, the existence of a PMF 
in our universe opens the opportunity to place indirect 
constraints on exotic particle physics models where the 
new physics couples to electromagnetism. In this paper 
we will investigate the consequences of a PMF for mod¬ 
els that contain magnetic monopoles, axions, and Dirac 
neutrinos with a magnetic moment. 

Blazars that emit TeV gamma rays are expected 
to produce an electromagnetic cascade of lower energy 
gamma rays due to electron-positron pair production and 
the subsequent inverse Compton up-scattering of cosmic 
microwave background (CMB) photons [10-14]. In the 
presence of an intergalactic magnetic held, electrons and 
positrons directed toward the Earth can be dehected off 
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of the line of sight, and those that are directed away from 
the Earth can be dehected back toward the line of sight. 
As a result the point source hux is depleted in the GeV 
band, and the blazar acquires a halo of GeV gamma rays. 
The non-observation of these GeV gamma rays was used 
to place a lower bound on the magnetic held strength at 
the level oi B> 10“^® G [15-17]. This bound depends on 
modeling of the blazar hux stability and also the plasma 
instabilities during propagation, and it may weaken sub¬ 
stantially depending on these assumptions [18-21]. The 
search for the GeV halo extended emission has been ongo¬ 
ing [22-26]. Most recently, Chen et al. [27] have found ev¬ 
idence for the halo in a stacked analysis of known blazars 
at ^ 1 GeV energies and interpret it to be due to a held 
with strength B ~ 10“^’’ — 10“^^ G. For reference, mea¬ 
surements of the cosmic microwave background place an 
upper bound on the magnetic held strength at the level 
of B < 10-9 G [28]. 

There are theoretical motivations for considering the 
possibility that the PMF is helical, i.e. there is an ex¬ 
cess of power in either right- or left-circular polarization 
modes. Helical magnetic helds emerge in many mod¬ 
els of magnetogenesis [5-7, 29] , and helicity conservation 
dramatically impacts the evolution of the PMF, aiding 
in its survival and growth [30]. Recently, Tashiro et al. 
[31, 32] analyzed the diffuse gamma ray sky at 10-60 GeV 
energies to look for the parity-violating signature [33, 34] 
of a helical magnetic held. They hnd evidence for an 
intergalactic magnetic held with strength B ~ lO”^"* G 
on coherence scales Xb ^ 10 Mpc and with left-handed 
helicity. Although the results of Refs. [31, 32] and [27] 
appear inconsistent, it may be possible to reconcile them 
by noting that the weak bending approximation breaks 
down for B ~ 10“^"^ G for gamma rays at ^ 1 GeV ener¬ 
gies [35]. 
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Motivated by these recent results, we investigate if the 
existence of a helical PMF can be used to constrain other 
particle physics ideas in a cosmological setting. Our anal¬ 
ysis is sufficiently general that our results will remain rel¬ 
evant even if the gamma ray observation results should 
change or go away, e.g. with more data. However, for the 
purpose of numerical estimates we will use B ^ 10“^^ G 
and As ~ lOMpc as the fiducial field strength and co¬ 
herence length scale, and we will take the magnetic field 
to have maximal (left-handed) helicity. 

In Sec. II we consider the interaction of a hypothetical 
abundance of cosmic magnetic monopoles and a PMF. 
The magnetic held does work on the monopoles, and 
its held strength is thereby depleted. The constraints 
we obtain in this way are generally weaker than exist¬ 
ing bounds but are obtained under a different set of as¬ 
sumptions. These results are summarized in Fig. 1. We 
also discuss how heavy magnetic monopoles can lead to 
anomalous scaling of the energy density in the PMF. As 
observations of the PMF improve, they can be sensitive 
to the anomalous scaling and thus become a tool for fur¬ 
ther constraining magnetic monopoles. 

In Sec. Ill we consider the interaction of an axion [(f) 
with the PMF through the coupling gaj(pFF. In this 
analysis we include the cosmological plasma, and thus 
we study the equations of magnetohydrodynamics cou¬ 
pled to an axion. Although we find that the axion has 
a negligible effect on the spectrum and evolution of the 
PMF, it is interesting to note that this conclusion is not 
sensitive to the assumed scale of Peccei-Quinn symmetry 
breaking, fa, as long as ga-y oc l//a- In turn, the PMF 
leaves the evolution of the axion condensate largely un¬ 
affected. In principle the PMF damps the axion oscilla¬ 
tions and the helicity of the PMF shifts the equilibrium 
point, but these effects are quadratic in the already-small 
magnetic held strength. 

In Sec. IV we consider the interaction of the neutrino 
magnetic dipole moment, with the PMF. Enqvist et 
al. [36, 37] have shown that this interaction induces a 
spin-hip transition, which cannot be in equilibrium in the 
early universe without running afoul of constraints on the 
number of relativistic neutrino species. Using their result 
with our hducial value of H ^ G, we evaluate an 

upper bound on /Xj^, which is shown in Fig. 2. 

We work in the GGS system with h = c = I. The unit 
of electric charge is e = y/a ~ 0.085 with a ~ 1/137 the 
hne structure constant, and the unit of magnetic charge 
is era = l/2e ~ 5.9. The magnetic held is measured 
in Gauss, and IG ~ 6.93 x 10“^°GeV^. The reduced 
Planck mass is denoted by Mp ~ 2.4 x 10^® GeV. The 
metric signature is (-1-), and the antisymmetric ten¬ 

sor normalization is = _|_x. 


II. MAGNETIC MONOPOLES 

A conservative cosmological bound on the energy den¬ 
sity of magnetic monopoles is flm = Pm/Pa < 0.3 


where pra is the energy density in monopoles and pcr — 
10“29gmc^/ cm® is the critical cosmological energy den¬ 
sity. The number density of nonrelativistic monopoles is 
nm = Pm/fn with m the monopole mass, and the cosmo¬ 
logical bound implies 

( 1 ) 

The bound grows weaker for lighter monopoles since they 
contribute less to the energy density for the same number 
density. 

The existence of the galactic magnetic held leads to 
another indirect bound. Magnetic monopoles tend to de¬ 
plete a magnetic held in the same way that free elec¬ 
trons short out a conductor. The survival of the micro- 
Gauss galactic magnetic held implies an upper bound on 
the directed hux, F, of magnetic charge onto the Milky 
Way. Requiring that the time scale for B-held depletion 
is longer than the dynamo time scale of B-held regenera¬ 
tion (rdyn — 10® yr), leads to the so called Parker bound 
[38] 


< 0.9 X 10 cm ^ sec ^ sr ^ . (2) 


Assuming that monopoles have unit charge and travel 
with velocity v, the Parker bound can be expressed as an 
upper bound on the monopole number density: 


(47r sr) F 


<(4x 10 2® cm ')(y^) 


-1 


( 3 ) 


Just as the Parker bound is predicated on the existence 
of a magnetic held in the Milky Way, we expect that 
a similar bound can be inferred from the existence of a 
primordial magnetic held in the early universe. 

We study a gas of monopoles and antimonopoles im¬ 
mersed in a magnetic held that permeates the cosmo¬ 
logical medium. The monopoles have mass m, magnetic 
charge Cm, and they are homogeneously distributed with 
number density nm{t). The magnetic held B(t,x) in¬ 
duces a Lorentz force of 


Fb = e^B (4) 

on a monopole at (t,x), and it begins to drift along the 
held line with a velocity v. The held does work by push¬ 
ing the monopole, and in this way the monopole extracts 
energy from the magnetic held at a rate 

Ub = Fb • V. (5) 

To solve for the evolution of the magnetic held strength 
we must know the monopole velocity. Prior to electron- 
positron annihilation, the monopole’s velocity is re¬ 
stricted by elastic scattering with the cosmological 
medium, but afterward it can free stream. We will con¬ 
sider each of these cases in turn. 
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A. Friction Dominated Regime 

In the epoch prior to e“'“e“ annihilation the cosmolog¬ 
ical medium was dense with electromagnetically charged 
particles. In this regime, monopoles interact with the 
medium through elastic scattering such ds M + 

M + with M the monopole. It is safe to assume that 
the monopole’s rest mass is much larger than the kinetic 
energy of particles in the plasma. This allows us to char¬ 
acterize the effective interaction with a drag force that 
takes the form [60] 


^drag — /dragV ■ (6) 

At the time of interest, the scatterers are relativistic with 
energy comparable to the temperature T of the plasma, 
and they are in thermal equilibrium with number density 
n ^ T^. For such a system the drag coefficient takes the 
form [39, 40] 


/drag (7) 

where gem{t) is the number of relativistic, charged de¬ 
grees of freedom in thermal equilibrium at time t and 
/3 is an 0(1) number related to the spin character and 
charge of the scatterers. We will drop /3 from this point 
onward since it is parametrically redundant with gem- 
Also note that gem and T depend on time, but this can 
be ignored on short time scales as compared to the Hub¬ 
ble time scale. 

The monopole’s equation of motion is 


The magnetic field’s response to this current is given 
by the magnetic analog of Ampere’s law, 

B = -47r Jm = -47remnmV . (11) 

where n,^ is the number density of monopoles (assumed 
equal to the number density of antimonopoles). Here 
we have used E = 0 since electric fields are screened 
due to the high electrical conductivity of the cosmologi¬ 
cal medium. For typical parameters, the inter-monopole 
spacing is small compared to the correlation length of the 
magnetic field, Um ^ Ab, and we can interpret and 
B as coarse grained quantities on this length scale. Then 
we insert v = ■Wterm from Eq. (9) into Eq. (11) to get the 
solution, 

B(t) = B(tj) e^^^^^AZ-r-dacay (']^2) 

where the decay time scale of the magnetic field is given 
by 


^decay — 




(13) 


In obtaining this solution, we have assumed v = Vterm 
which is justified if the monopoles reach terminal veloc¬ 
ity much more quickly than the decay time scale, i.e. 
Tdecay Tterm- As we will See below, this condition is 
satisfied for the range of parameters of interest to us. 

To ensure survival of the magnetic field, we require 
that T(jecay is much larger than the Hubble time at tem¬ 
perature T, 


mV = em'B - /dragV. 


( 8 ) 


We assume that the distance traveled by the monopole 
is small compared to the correlation length As of the 
magnetic field, and we treat B as uniform. Eq. (8) im¬ 
mediately gives the terminal velocity of the monopoles, 


^term 


,B 




J^2 


which is achieved on a time scale 


Aerm 


m 

p2p2 „ ^ ■ 


( 9 ) 


( 10 ) 


Comparing with the Hubble time tn ^ Mp/T'^ (radia¬ 
tion era) we have Tterm ^ t-H provided that m < gemMp, 
and thus the cosmological expansion is negligible. 

At the present cosmic epoch, the photon temperature 
is ~ 10“'^ eV and B ~ 10“^"^ G. Assuming B oc T^, we 
get B ~ 10® G at T ~ MeV when gem — 10. These esti¬ 
mates give Vterm — 10“®, which validates our uses of the 
non-relativistic equation of motion. As a consequence, 
the distance traveled by a monopole during Tterm is quite 
small, dterm < 10“®Tterm- We shall assume that the cor¬ 
relation length of the magnetic field is larger, Ap > dtermj 
thus justifying our treatment of the magnetic field as be¬ 
ing uniform. 


" 2H ■ 


(14) 


with H the Hubble parameter and g* the effective num¬ 
ber of relativistic degrees of freedom. Substitution of 
Eq. (13) into r^ecay > tp now leads to a constraint on 
the number density of monopoles. 


nm < 


gemg*T^ 
67re^ Mp 


(15) 


when the universe had temperature T. 

The strongest bound is obtained when T is smallest. 
Since our calculation assumes that monopoles scatter 
on relativistic, charged particles with a thermal abun¬ 
dance, the last time at which this is possible is the 
epoch of e+e” annihilation. At this time Tann ~ 1 MeV, 
gem ~ g* ~ g*s — 10.75. To translate this into a 
bound on the monopole number density today, denoted 
by no, we multiply by (oann/oo)® = (g*S,o7o/g^S'T’ami) 
with g-tSfi — 3.91 and Tq ~ 2.3 x 10“^ eV the temper¬ 
ature of the microwave background photons today. We 
find an upper bound on the monopole number density 
today, 


no < 


^ ^m g*5,0g*gem '^ann^Q 
Mp 


~ 1 X 10"2°cm-®. (16) 
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FIG. 1: A summary of upper bounds on the magnetic 

monopole abundance from this work and the literature. 
Black: the requirement of survival of the primordial magnetic 
field (“primordial Parker bound”), derived here in Eq. (16). 
Red: the cosmological abundance bound in Eq. (1). Blue: 
direct search constraints [42] (see the magnetic monopole re¬ 
view). Green: the requirement of survival of the Galactic 
magnetic field (Parker bound), given by Ref. [43]. Orange: 
the requirement of survival of the Galactic seed field (“ex¬ 
tended Parker bound”), given by Ref. [41]. We take v ~ 10“^ 
and assume that monopoles are unclustered, fc — 1. If the 
monopoles are clustered then the Parker bound, extended 
Parker bound, and direct search limits move down by a factor 
of/, ~ 10®. 


If this bound is not satisfied then any primordial mag¬ 
netic field would have been exponentially depleted by the 
time of electron-positron annihilation. Due to the close 
connection with the Parker bound for survival of galactic 
magnetic fields, we will refer to Eq. (16) as the “primor¬ 
dial Parker bound.” 

In Fig. 1 we compare the primordial Parker bound 
in Eq. (16) with other constraints derived previously in 
the literature. Since these constraints are typically ex¬ 
pressed as a bound on the monopole flux, we translate 
into a bound on the number density using (47rsr)J^ « 
fcTiv. Here v is the average monopole velocity and 
fc = ?^gaiaxy/?^cosmo. is the enhancement factor that ac¬ 
counts for clustering of monopoles in the galaxy. For 
clustered monopoles fc ~ 10®, but otherwise /c 1. 
For the extended Parker bound calculation of Ref. [41], 
we take Bseed = 10“^^ G. For the direct search con¬ 
straint we show a relatively conservative and robust limit 
of < 10“^® cm“^ sec“^ sr“^, but stronger constraints 
are available for specific monopole parameters [42] . From 
the figure, one can see that the primordial Parker bound 


becomes stronger than the cosmological bound for light 
monopoles, m < 5 x 10^® GeV, but it always remains 
weaker than the direct search constraints. 


B. Free Streaming Regime 


After cosmological electron-positron annihilation the 
number density of these scatterers decreases by a factor 
of ^ 10“®®. The monopoles experience very little drag 
force, and they can be accelerated freely by the magnetic 
field. For a uniform and static magnetic held, the solution 
of Eq. (8) with /drag = 0 is simply v = emBt/m, or for 
an inhomogenous held with domains of size \b we hnd 


v{t) 


m 



(17) 


if the motion is diffusive. The monopole becomes rela¬ 
tivistic when v{tcei) 1 and comparing this time with 
the present age of the universe gives 


trei ^ / m y to ^ 
to \emBtoJ Xb \MpJ Xb 


With Xb ~ Mpc and to ^ lOGpc, we hnd that 
monopoles are relativistic today if to < 10® GeV, and 
they are non-relativistic otherwise. 

The above estimate ignores backreaction of the 
monopoles on the PMF. To check for consistency, we 
compare the kinetic energy in monopoles pkin to the en¬ 
ergy density available in the PMF pB = B^lSn. For rel¬ 
ativistic monopoles we should have pb > Pkin ^ 'mrim, 
and this provides an upper bound on the number density 
of monopoles for which the velocity estimate in Eq. (17) 
can be expected to hold. Taking B ^ 10“®"® G we hnd 


Um < 10 


-35 


10® GeV \ 

TO J 


cm 


-3 


(19) 


Since are interested in much larger number densities, as 
indicated by the bound in Eq. (16), we cannot use the 
velocity relation in Eq. (17), but instead the monopole 
and magnetic held equations will need to be evolved si¬ 
multaneously. 

Without friction to provide a means of energy dissi¬ 
pation, the monopoles cannot deplete the magnetic held 
strength. Instead there is a conservative exchange of en¬ 
ergy between the magnetic held and the kinetic energy of 
the monopoles. This co-evolution can lead to an anoma¬ 
lous departure from the usual power law scaling behav¬ 
ior of the magnetic held energy density if the monopoles 
are non-relativistic. This can be seen from the following 
argument. In the absence of the monopole gas, the en¬ 
ergy density in the magnetic held redshifts like radiation 
Pb ~ (1 + -^)"^ where 2 ; is the cosmological redshift. Mean¬ 
while the kinetic energy stored in a gas of non-relativistic 
particles redshifts more quickly. We can write the kinetic 
energy density as pkin = (2 to) where p is the typi¬ 

cal momentum and rim is the monopole number density. 
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Since p ^ {1 + z) and Um ^ (1 + z)^, the kinetic energy 
density redshifts like pkin ^ (1 + z)^. If energy is trans¬ 
ferred quickly between the monopoles and magnetic field 
then we might expect ps Pkin. ~ (1 + z)®/^ where 9/2 
is the average of 4 and 5. This result is confirmed by the 
full calculation, which we will now present. 

Once again we consider a gas of monopoles coupled 
to a magnetic field, but we now include the effects of 
cosmological expansion. Let X^(t) be the world line of 
a monopole, and let U^(t) = dX^/dr be its 4-velocity. 
The monopole equation of motion, given previously by 
Eq. (8), is now replaced by 

= e^F\Ui^ (20) 

where T^^ is the Christoffel symbol. The magnetic analog 
of Ampere’s law in Eq. (11) is now replaced by 

= 47r4 (21) 

where Vq, is the covariant derivative. The magnetic cur¬ 
rent density arises from monopoles with velocity C/^ 
and antimonopoles with velocity t/^. It can be written 
as 


m 


dU^ 

dr 


I 

' rr 


a/Hiof, and the solution is 


B = 


Vm = 
where 


-9/4 


-3/4 


Jo(</>) Bi+yo(,/)) B2 

(a/ai)-i/4 


(p Ji((/) Vi+(j) Fi((/) V2 





a* 


(26a) 

(26b) 

(27) 


At late times p ^ 1, and all of the Bessel functions go 
to zero with an envelop ~ ~ Then the 

terms in square brackets do not scale with a and we find 
B ~ and Vm ~ One can check that the 

energy densities scale in the same way 


|B| 


Pb = 


oc a 


-9/2 


Stt 

^ I i2 —Q/2 

/^kin. — ^ ^ ; 


(28a) 

(28b) 


which confirms our earlier argument. 

During the matter era we have rj = rji + — 

« 2yaj, and the solution is 


Jm = e™ {rim - rifn U^) , (22) 

and it satisfies the conservation law 

V^^ = 0. (23) 

The spatial component of the 4-velocity is the comoving 
peculiar velocity (11^)* = Uy, and with an additional 
factor of a we form the physical peculiar velocity (v^)* = 
ally. Neglecting the electric field and spatial gradients, 
the system of equations can be put into the form 


dr){a^m) = —(o^B) 
m 

(24a) 

a<9^(a^B) = -^rem[o^nm) (av^) 

(24b) 

d-q (^U = 0 

(24c) 


where dp = dt/a = da/Ha^ is the conformal time coor¬ 
dinate. The third equation implies that the number of 
monopoles per comoving volume is conserved; no = a?nm 
is the number density of monopoles today. Then if not 
for the additional factor of a in the second equation, the 
solutions would simply be oscillatory with angular fre¬ 
quency 


B = 


-9/4 


E 

s=±l 


Isy/AQ’^y/ai-l 


B. 


v„, = I — 




s=±l 


(29a) 

(29b) 


For typical parameters we have ACo'^pflai 1, and the 
solution is oscillatory with a power law envelope. We find 
the same anomalous scaling as in the radiation era, c/. 
Eq. (28). If we remove the monopoles from the problem 
by sending nm,dj 0 then the would-be oscillatory fac¬ 
tors become a power law decay, and we regain the usual 
scaling B ~ and Vm ~ a~^. 

This anomalous scaling does not provide constraints on 
the monopole number density. However, it does affect the 
way that we translate constraints on the magnetic field in 
the early universe into the value of the magnetic field to¬ 
day. Measurements of the cosmic microwave background 
restrict the magnetic field energy density to be less than 
~ 10“® of the photon energy density at the time of re¬ 
combination [28] 

PBiZrec) < 10"^P-y(Zrec) ■ (30) 


. /47re^no 

This is just the usual formula for plasma frequency, but 
instead of electron charge, mass, and density, here we 
find the corresponding parameters for the monopole gas. 

To solve these equations we must relate a to p. Dur¬ 
ing the radiation era we have p = pi + {a — ay/Hiaf « 


Using the scaling relations, ps ~ (1 -I- z)®/® and p^ ^ 
(1 -|- z)4 this inequality implies that the magnetic field 
energy density today is bounded by 


PBfi < 10 ^P 7 ,o(l + Zrec) 


~ (3 X 10"^° G)^ 


1 + 


- 1/2 


1300 


(31) 
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where ~ 2-k'^Tq/ZQ is the CMB energy density to¬ 
day. If it were not for the anomalous redshifting, the 
constraint on the B-field strength would be weaker by a 
factor of (1 -I- — 6. 

The anomalous scaling of the field strength in Eq. (28) 
can become a tool in the future as measurements of the 
PMF improve. By measuring the magnetic field strength 
at different redshifts, say by using TeV blazars at differ¬ 
ent distances, we can directly probe the anomalous scal¬ 
ing, and hence obtain a new handle on the relic density 
of non-relativistic magnetic monopoles. 


III. AXIONS 


Consider an axion coupled to the electromagnetic 
field A^(x). The Lagrangian takes the form [44] 


IbTT^ 


(32) 


where F^'' = F^p is the dual field strength ten¬ 

sor, and = (p, j) is the electromagnetic current arising 
from the charged Standard Model fields. Our analysis 
is sufficiently general to apply to any axion or axion-like 
particle described by Eq. (32), but as a fiducial refer¬ 
ence point we will consider a QCD axion with Peccei- 
Quinn scale of fa — 10^° GeV, an axion mass of rua « 
^qcd/ fa — ImeV, and a photon-axion coupling constant 
gaj « a/(27r/a) ~ lO-^^GeV"^ 

The classical axion condensate obeys the field equation 


Uip + 


9ay 

Att 


E B 


(33) 


We seek to study the coevolution of the coupled ax¬ 
ion and electromagnetic fields. Eqs. (33) and (36) de¬ 
scribe a non-dissipative system. Dissipation is introduced 
as the electromagnetic field couples to charged particles 
in the cosmological medium, which opens an avenue for 
energy to be lost in the form of heat. This coupling 
is parametrized by the conductivity a which appears in 
Ohm’s law 


j = a(E + vxB) (37) 

where v(t,x) is the local velocity of the plasma. Prior 
to the epoch of e“'"e“ annihilation, free charge carriers 
were abundant and the cosmological medium had a high 
conductivity [45] 


cr « T/a (38) 

where a ~ 1/137 is the fine structure constant. Ohm’s 
law allows us to eliminate j and thereby reduce the sys¬ 
tem of equations in four unknowns {E, B,j,(p} to a set 
of equations describing only three unknowns: 

if - VV + rnlip = ^E • B (39a) 

47r 

B = - V X E (39b) 

E = V X B - gajfB - gajVip x E 

— dTTcrE — 47rcrv X B . (39c) 

In the MHD approximation (nonrelativistic flow) we can 
neglect the displacement current since it is negligible 
compared to the curl of the magnetic field, |E|/| V x B| ^ 
{v/cY ^ 1 [46]. Then Eq. (39c) becomes algebraic in E, 
and we can solve it to eliminate E from the remaining 
equations. Focusing now on a homogenous axion field, 
the system of equations reduces to 


where we have used F^^F^^’-' = —4E-B. The electromag¬ 
netic field evolves according to the constraint equation 

= 0 (34) 

and the modified field equation 

+ga^d^^F>^'^ = . (35) 

In terms of the electric and magnetic vector fields we have 

V • B = 0 (36a) 

1 /9R 

VxE-H-—= 0 (36b) 

c at 

V • E -I- ga-y^f ■ B = 47rp (36c) 

1 i9E 47r 

^ y-'B-- gayfB-gayV(pxE=—j. (36d) 


... 2 ^djB] 


-If + m^ip = • V X B 


B = V X (v X B) -I- 77d V^B gayVdf V x B 


(40a) 

(40b) 


where 

_ 1 a 

4707 47rT 


(41) 


is the magnetic diffusivity, assumed to be homogenous. 
Eqs. (40a) and (40b) together with the Navier-Stokes 
equation for the plasma velocity v are the final equations 
to be solved. We first solve Eq. (40b) to determine the 
effect of the axion on the magnetic field, and afterward 
we will consider the evolution of the axion according to 
Eq. (40a). 


A. Effect of Axion on Magnetic Field 


Note the presence of the additional terms arising from 
the spatio-temporal variation of the axion field. 


In order to solve Eq. (40b) for the B-field we must know 
the fluid velocity, which appears in the advection term. 
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V X (v X B). Since our primary interest is in the co¬ 
evolution of the magnetic field and the axion, not in the 
magnetohydrodynamics, we will neglect this term [61]. 
Since Eq. (40b) is linear in B, and we assume that tp 
is homogenous, we can solve the equation by first per¬ 
forming a Fourier transform. For a given mode k let 
(ei(k),e 2 (k),e 3 (k)) form a right-handed, orthonormal 
triad of unit vectors with e 3 (k) = k/jkj. It is convenient 
to introduce the right- and left-circular polarization vec¬ 
tors by 


„±n ^ ei(k)±ie2(k) 

e (k)-. 


(42) 


Note that ik x e^(k) = ±|k|e^(k). The mode decompo¬ 
sition is given by 


/ lj3 h 

^e**^-^6,(t,|k|)e^(k). (43) 


With this replacement, Eq. (40b) becomes 

h± (t, k) = -r]dk'^b± {t, k) ± kipb± {t, k) (44) 

where we have written k = jkj. The last term of this 
equation, essentially the chiral-magnetic effect [47], has 
been studied previously in the context of axions [48] and 
cosmology [49, 50]. The solution is 

b±{t,k) = b±iU,k) 

where we have defined the wavenumber 

2 

and A(p(t) = <p(k) — <p[ti) is the change in the axion field. 
The prefactor in Eq. (45) is the initial spectrum of the 
magnetic field which will depend on the PMF genera¬ 
tion mechanism. The first exponential is the usual diffu¬ 
sive decay term, which exponentially suppresses modes 
on a length scale shorter than = \J[t — ti)r]d — 

10“^ \/tjT. The second exponential only kicks in at small 
length scales where k > k^x- Then it leads to a suppres¬ 
sion of one polarization mode and enhancement of the 
other, depending on the sign of A<p(t). 

The value of A(p{t) depends on the solution for the 
axion field as well as the initial time ti. We expect that 
the misalignment mechanism sets the initial condition 
ip{ti) ^ ±fa- The subsequent evolution is determined 
by solving Eq. (40a), which we will turn to in the next 
section. For the moment we will assume that the axion 
evolution is not significantly affected by the presence of 
the magnetic field, and the solution is the standard one: 
the axion remains “frozen” at ip(ti) until the time of the 
QCD phase transition when it begins to oscillate around 
ip = 0 with angular frequency uj = rua [51]. Then we can 
approximate 

A„(() =» 


t < t 
t > t 


QCD 


(47) 


QCD 


(45) 

(46) 


where s = sign[(/3(ti)]. Using this approximation we can 
estimate k^^- Prior to the QCD phase transition, Aip is 
small and fcax is large, meaning that none of the modes 
receive the enhancement or suppression from the axion 
coupling. This is reasonable since the axion is deriva¬ 
tively coupled, and as long as it is stationary there will 
be no effect on the magnetic field. After the QCD tran¬ 
sition, we can estimate 


k 


ax 


Aira 


Qa'yfa 


AirT 


(48) 


using Eqs. (38) and (41). Note that this result is insen¬ 
sitive to the Peccei-Quinn scale, and as long as ga-y = 
aj{2Trfa) we have k~J^ « ~ 10“®T“^. 

The solution in Eq. (45) can also be written as 


b±{t,k) = b±{U,k) e-(fcT^)Ut-t.Ud (49) 


where 


K{t) 


1 

^ax(^ 


ga^Aipjt) 
2{t-U) ■ 


(50) 


This representation of the solution is convenient, because 
all the spectral information is contained in the second 
factor. One of the helicity modes has a G aussian spec- 
trum peaked at \K{t)\ > 0 with width 
and the other helicity mode peaks at fc = 0. Estimating 
fcax as above, we find that the associated length scale of 
the spectral peak corresponds to K~^ « k^ytrid ~ 600t, 
which is larger than the scale of the cosmological horizon 
dn ^ t. 

It appears that the presence of an axion condensate 
coupled to electromagnetism has a negligible impact on 
the evolution of a primordial magnetic field, unless there 
are situations in which Aip can be much larger than fa- 

We note that our analysis ignores the possibility of 
turbulence in the primordial plasma. It would be of in¬ 
terest to include both turbulence and the axion coupling 
in future studies. 


B. Effect of Magnetic Field on Axion 


Next we will investigate the effect of a background 
magnetic field on the axion condensate. We have seen 
that the magnetic field is approximately unmodified on 
length scales larger than the diffusion length, fcjjg ^ 
{cf., Eq. (45)). In this regime Eq. (40a) can be 
rewritten as 


where 


p -1- 2 - \- — = % 

(51) 

Aiecay Ti 

Stt 

(52) 


I 

Ta= - , 

(53) 

rria 


n- ^7'^‘^(B-V X B), 

47r 

(54) 














and the angled brackets (•) denote spatial averaging. The 
axion condensate evolves like a damped and driven har¬ 
monic oscillator, where the damping and driving forces 
are induced by the magnetic field background. As we 
discuss below, it is interesting that the driving force is 
associated with the helicity of the magnetic field. 

The magnetic-induced damping of axion oscillations 
is parametrized by the time scale Tdecay To determine 
when this damping will be relevant for the evolution of 
the axion, we compare it with the cosmological time 
scale, given by Eq. (14). To express (|Bp) = in 
terms of the magnetic field strength today, Bq, we use 
B = Bo{ao/a)'^ ~ 10Bo{T/T q)"^ where the factor of 10 is 
related to the number of relativistic degrees of freedom 
in the early universe and today. Then the ratio is found 
to be 


We are led to conclude that for realistic parameters, the 
helical PMF does not significantly impact the evolution of 
the axion condensate prior to the QCD phase transition. 

After the QCD phase transition, the axion mass 
reaches its asymptotic value, and the source term dis¬ 
places the minimum of the axion potential from ip = 0 
to (pmin = = H/m^. In terms of the angular coordi¬ 

nate: 


fa ^ 

10-4" 


25 agq-yBlT^ 
127r2 famlXBfiT^ 


Bp \ ^ / Xb,o 
10-14 cy VlOMpc 


-1 


( /a \ ( 9ai 

VlOiOGeVy 'v 10-13 GeV-i 


(57) 


200 MeV J 



'^decay ^ IOtT Tq 

^ 75 agl^MpB^T ^ > 

^^16(10-13 GeV-4)2 (lOio GeV) (10 -i4G)2 
~ 10 - ^ .. ^ ~ . 

i/ay ^ ^P 

This estimate suggests that the magnetic-induced decay 
of the axion field is negligible for a typical Peccei-Quinn 
scale and B-field strength. If the B-field strength today 
were as large as Bp ~ 10-® G and the Peccei-Quinn scale 
was as low as fa TeV, then Tdecay would be compara¬ 
ble to the Hubble time at T « /„. As the temperature 
decreases, the magnetic-induced decay becomes less rel¬ 
evant. 

It is interesting that the magnetic field also induces a 
driving force, parametrized by B. The pseudoscalar B 
is related to the helicity of the magnetic field. This is 
perhaps more evident from the initial form of the axion 
field equation, Eq. (33), where E • B is equal to the rate 
of change of the helicity density —(l/2)d(A • li)/dt plus 
a divergence, which vanishes upon spatial averaging. If 
the power in the magnetic field is localized on a particular 
length scale Ab we can estimate (B •VxB)~BVAb~ 
300{Bq/Xbo){T/T p)^ where we used Xb ~ 3A_b,o(To/T). 

Prior to the QCD phase transition we can neglect the 
mass and drag terms in Eq. (51), and the solution is 
simply p = Bt^ 12. Since the axion is massless, there is 
no restorative potential, and the helical magnetic field 
leads to an unbounded growth of the axion condensate. 
Although this analysis neglects the Hubble drag, we can 
estimate the maximum field excursion in one one Hubble 
time to be Bt'jj/2. Comparing with the Peccei-Quinn 
breaking scale, the corresponding angular excursion is 

A g ^ ^ 75 agq^MlBl 

fa 167r2 faglTpXB,P 

^10-33 

V10-14 g; VioMpcy 


The temperature dependence enters through B ^ T^, 
Xb ~ 1/T, and rjd ^ 1/T, and the fractional shift is 
largest at high temperature where B is large and Xb 
is small. Immediately after the QCD phase transition, 
Tqcd ~ 200 MeV, the fractional shift is already extremely 
small. Moreover if the PMF is not helical then B = 0 
and there is no shift in the axion potential. 

It is interesting that the estimate of Eq. (57) is insensi¬ 
tive to the Peccei-Quinn scale; as long as = otl^T^fa) 
and TOa = Xff^calfa we have ga-y/faml = Q;/(27rA4pj,). 
Then the primarily challenge toward obtaining a large 
effect is the smallness of the magnetic field strength. Al¬ 
though unrelated to primordial magnetic fields, which is 
the motivation for this work, it would be interesting to 
study the axion condensate in an astrophysical system 
where the magnetic field is both helical and strong. For 
instance, the field strength in a magnetar can grow as 
large as B ~ lO^^G and the magnetic field in some as¬ 
trophysical jets is known to be helical [52]. 


C. Axion-Photon Interconversion 

Until this point we have focused our attention on the 
interplay between the axion condensate and the primor¬ 
dial magnetic field, and we now turn our attention to the 
quanta of these fields. In the presence of a background 
magnetic field, the interaction in Eq. (32) yields a mix¬ 
ing between axion particles and photons. Typically the 
conversion is inefficient, but in the presence of a plasma 
the photon acquires an effective mass, and the conver¬ 
sion probability experiences a resonance when m.y = rua 
[53]. In the cosmological context, the conversion of pho¬ 
tons into axions may lead to a dimming of the cosmic 
microwave background across frequencies. Then mea¬ 
surements of the spectrum of the GMB can be used to 
place constraints on the axion-photon coupling and the 
magnetic field strength. 

Bounds were obtained from the COBE / FIRAS mea¬ 
surement of the GMB spectrum in Ref. [54], and re¬ 
cently a second group [55] has extended the calculation 
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to include forecasts for next-generation CMB telescopes, 
namely PIXIE and PRISM. The latter references finds 
an upper bound on the product of the axion-photon cou¬ 
pling and the r.m.s. magnetic field strength today: 

gajBo < GeV^ nG (COBE-FIRAS data) 

gajBo < 10"^® GeV^ nG (PIXIE/PRISM forecast) 

(58) 

for a light axion nia < 10“eV. For larger axions masses 
the bound weakens. Using our fiducial value for the mag¬ 
netic field Bq ~ 10“^"* G, we can write 

gay < 10-9 GeV-i (^ 0 ^^) (GOBE-FIRAS) 

gay < 10-“ GeV-i (PIXIE/PRISM). 

(59) 

These bounds are comparable to the direct search limits 
from the GAST helioscope [56], 

gay < 8.8 X 10““ GeV“^ (60) 

for TTia < 0.02 eV. 

IV. DIRAC NEUTRINOS 

While the neutrinos are known to be massive particles, 
the nature of their mass remains a mystery. If neutrinos 
are Dirac particles then the theory contains four light 
states per generation: an active neutrino vl, an active 
antineutrino. Dr, a sterile neutrino vr, and a sterile an¬ 
tineutrino vr. The active states interact through the 
weak force, and this allows them to come into thermal 
equilibrium in the early universe. The sterile states, on 
the other hand, interact only via the Yukawa interaction 
with the Higgs boson, and because of the smallness of 
the Yukawa coupling ~ rua/v ~ these states 

are not expected to be populated. 

This story is modified if a strong magnetic field per¬ 
meated the early universe. The nonzero neutrino mass 
implies that the neutrino will also have a nonzero mag¬ 
netic moment /i^. From Standard Model physics alone 
one expects [57, 58] 

(61) 

where gs = ej^ma ~ 83.6 GeV“^ is the Bohr magne¬ 
ton, but new physics can increase this value appreciably. 
The magnetic field couples to and induces the spin- 
flip transitions vr vr and Dr Dr, which can be 
viewed as the absorption or emission of a photon. If the 
spin-flip occurs rapidly in the early universe, the sterile 
states would be populated, and the effective number of 
relativistic neutrino species would double from W = 3 to 


6. However, this is not consistent with measured abun¬ 
dances of the light elements, which imply N^, « 3 at 
the time of nucleosynthesis [59]. We must therefore re¬ 
quire that the spin-flip transition goes out of equilibrium 
prior to the QCD epoch, Tqcd — 200 MeV, so that the 
subsequent entropy injection at the QGD phase transi¬ 
tion can suppress the relative abundance of sterile states 
to acceptable levels [51]. This translates into an upper 
bound on the neutrino magnetic moment and magnetic 
field strength, which was originally discussed by Enqvist 
et al. [36, 37]. 

In the rest of this section we apply the results of 
Ref. [37]. The spin-flip transition occurs with a rate 

Tr^r = (62) 

where is the average conversion probability 

and is the total weak scattering rate. The active 
neutrinos scatter via the weak interaction which leads 
to ~ 30 G|.Tqpj 3 at the QGD epoch. The conver¬ 
sion probability depends on the magnetic moment and 
field strength as cx g^B'^, since the interaction 

Hamiltonian is i7int = —Mi/ • B. The coefficient takes 
different values depending on the relative scale of the 
magnetic field domains Xr and the weak collision length 
Lw « (r(^*)-^. At the QGD epoch Lw — 1-6 x IQ-^ cm, 
which corresponds to a length scale of Lwfi — 3 x 10^° cm 
today. It is safe to assume that the magnetic field of in¬ 
terest is much larger than this length scale, and therefore 
As 3> Lw- To ensure that the spin-flip transition is out 
of equilibrium one must impose < H with H the 

Hubble parameter. This inequality resolves to the bound 
(see Eq. (37) of Ref. [37]) 

< (3.5 X 10Vs G) . (63) 

To express this inequality in terms of the B-field strength 
and correlation length today, we use B ~ 6i?o(TQCD/7o)9 
and Lw/Xb = Lw,o/X b,o- This leads to an upper bound 
on the neutrino magnetic moment: 

(64) 

If this bound is not satisfied, the sterile neutrino states 
will still be thermalized with the active neutrino states at 
the time of BBN leading to W ~ 6, which is inconsistent 
with the data. If the neutrinos are Majorana particles, 
then the sterile states are much heavier, and this bound 
does not apply. 

The bound in Eq. (64) is represented graphically in 
Fig. 2. For comparison we show the SM prediction from 
Eq. (61) and the direct search limits. The strongest labo¬ 
ratory constraints arise from elastic v — e scattering. The 
limits are flavor-dependent, but they are typically at the 
level of [42] 

Dv ^ 10“^9/rB (direct). (65) 








10 



Field Strength Today: Bq [ Gauss ] 


FIG. 2: The requirement that spin-flip transitions are out of 
equilibrium at the QCD epoch leads to an upper bound the 
neutrino magnetic moment given by Eq. (64). 


From the figure we see that the indirect early universe 
constraint is significantly stronger than the direct con¬ 
straint for B > 10“^® G. This provides the exciting op¬ 
portunity to constrain extensions of the SM that predict 
an enhancement to the magnetic moment of Dirac neu¬ 
trinos. 


V. SUMMARY 

Growing evidence for the existence of an intergalactic 
magnetic field has motivated us to consider the effects of 
a primordial magnetic field on models of exotic particle 
physics in the early universe. We have focused our study 
on magnetic monopoles, axions, and Dirac neutrinos with 
a magnetic moment. We summarize our results here. 

In the context of a universe containing relic magnetic 
monopoles, we have derived a “primordial Parker bound” 
by requiring the survival of a primordial magnetic field 
until the time of electron-positron annihilation. The 
bound, which appears in Eq. (16), gives an upper limit on 
the cosmological monopole number density today: uq < 
1 X 10“^° cm“®. This translates into an npper bound on 
the monopole flux in the Milky Way; if the monopoles are 
unclustered then <3x10 cm ^ sec ^ sr ^(n/10 ®), 
and if they are clustered the bound weakens by a factor of 
~ 10®. In Fig. 1 we compare the primordial Parker bound 
with other constraints on relic monopoles. If the primor¬ 
dial magnetic field is not generated prior to T ~ MeV, 
then this bound does not apply. 


After e+e” annihilation the monopoles are able to free 
stream, and they evolve along with the magnetic field as 
described by the system of equations in Eq. (24). The 
solution is an analog of the familiar plasma oscillations 
(“Langmuir oscillations”) seen in an electron-ion plasma. 
In the regime where the plasma oscillations are fast com¬ 
pared to the cosmological expansion, the coupling of the 
monopoles to the magnetic field leads to an anomalous 
scaling with redshift such that B ^ Vm a”®/"*, 

and ps ~ pkin. ~ The behavior of the coupled 

system is effectively the average of the usual scalings for 
radiation ps a~‘^ and the kinetic energy of a non- 
relativistic gas pkin. ~ If the strength of the inter¬ 
galactic magnetic field could be measured over a range of 
redshifts, this would allow for a direct test of the anoma¬ 
lous scaling, and thereby probe relic magnetic monopoles. 

We have also studied the effect of a primordial mag¬ 
netic field on the evolution of an axion condensate in the 
early universe. We obtain an exact solution to the MHD 
equations for the magnetic field in the limit where the ad- 
vection term is negligible and the axion is homogenous. 
After Peccei-Quinn breaking but prior to the QGD phase 
transition, the axion field is frozen, because its mass is 
smaller than the Hubble scale, and since the axion is 
derivatively coupled, this leads to no effect on the mag¬ 
netic field. Below the QGD scale the axion field begins 
to oscillate, and the spectrum of the magnetic field is 
distorted as in Eq. (45). One helicity mode of the mag¬ 
netic field is enhanced while the other is suppressed; this 
GP-violation is a consequence of the axion’s pseudoscalar 
nature. However, the spectral shape of the magnetic field 
is only affected on extremely large length scales, as given 
by Eq. (50), except in situations where there can be sig¬ 
nificant axion evolution prior to the QGD epoch. 

We next study the evolution of the homogenous ax¬ 
ion condensate in the presence of a background mag¬ 
netic field. The axion behaves as a damped and driven 
harmonic oscillator, as seen from its equation of mo¬ 
tion Eq. (51). The damping time scale depends on the 
strength of the magnetic field and the photon-axion cou¬ 
pling. Eor typical parameters it is generally larger than 
the cosmological time scale, and therefore irrelevant for 
the evolution of the axion. It is interesting that the driv¬ 
ing force (source term B) is only operative when the mag¬ 
netic field has a helicity. This can be seen directly from 
the interaction £ 9 (pFF where FF ^ E • B ~ is re¬ 
lated to the rate of change of magnetic helicity h = A • B. 
Prior to the QCD phase transition when the axion was 
effectively massless, the axion field equation reduces to 
(p = H. In principle a very strong magnetic field could 
cause the axion to grow without bound as ip{t) — Ht^/2, 
by drawing energy from the magnetic field. For typical 
parameters, however, this growth occurs on a time scale 
that is much longer than the cosmological time. It may 
still be the case that helical magnetic fields occurring in 
astrophysical environments are strong enough to lead to 
observable signatures. 

We have also considered the resonant conversion of 
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CMB photons into axions, which leads to a distortion of 
the CMB blackbody spectrum [53-55] . Using constraints 
on spectral distortions from current and anticipated fu¬ 
ture CMB telescopes, Ref. [55] obtained an upper bound 
on the axion-photon coupling. For our fiducial magnetic 
field strength this translates into ga-y ^ 10“® GeV“^ with 
current data, and a forecast of ga-^ < 10“^^ GeV~^ for 
experiments presently under discussion (see Eq. (59)). 

Finally we turn to the effect of the primordial mag¬ 
netic field on Dirac neutrinos, which carry a magnetic 
moment. In the presence of a magnetic field, left-handed 
neutrinos can be converted into right-handed neutrinos. 
If this spin-flip process is in equilibrium in the early uni¬ 
verse, the right-handed states would be populated, and 
the effective number of relativistic neutrino species would 
double from 3 to 6, which is inconsistent with observa¬ 
tions. Requiring that this process is out of equilibrium at 
the time of the QCD phase transition leads to an upper 
bound on the neutrino magnetic moment and magnetic 


field strength. Drawing on the work of Ref. [37], we find 
the limit in Eq. (64), which implies < 3 x 
for our fiducial magnetic field parameters Bq = 10“ G 
and Xb — lOMpc. As seen in Fig. 2, this bound is signif¬ 
icantly stronger than the direct search limits over most 
of the parameter space. 
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